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Abstract Dust is an important constituent in cometary comae; its analysis is 
one of the major objectives of ESA's Rosetta mission to comet 67P/Churyu- 
mov-Gerasimcnko (C-G). Several instruments aboard Rosetta are dedicated 
to studying various aspects of dust in the cometary coma, all of which re- 
quire a certain level of exposure to dust to achieve their goals. At the same 
time, impacts of dust particles can constitute a hazard to the spacecraft. To 
conciliate the demands of dust collection instruments and spacecraft safety, 
it is desirable to assess the dust environment in the coma even before the 
arrival of Rosetta. We describe the present status of modelling the dust coma 
of 67P/C-G and predict the speed and flux of dust in the coma, the dust 
fluence on a spacecraft along sample trajectories, and the radiation environ- 
ment in the coma. The model will need to be refined when more details of 
the coma are revealed by observations. An overview of astronomical observa- 
tions of 67P/C-G is given and model parameters are derived from these data 
where possible. For quantities not yet measured for 67P/C-G, we use values 
obtained for other comets. One of the most important and most controversial 
parameters is the dust mass distribution. We summarise the mass distribu- 
tion functions derived from the in-situ measurements at comet IP/Halley in 
1986. For 67P/C-G, constraining the mass distribution is currently only pos- 
sible by the analysis of astronomical images. We find that the results from 
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such analyses arc at present rather heterogeneous, and we identify a need to 
find a model that is reconcilable with all available observations. 



1 Introduction 

The inner comae of comets are the most dust rich environments in the solar 
system. Almost everything we see from a comet with the naked eye is dust. 
Both the coma and the tail are seen as sunlight scattered by /J,m-sized dust. 
Meteors, especially in meteor streams, are caused by mm- to cm-sized parti- 
cles that originate to a large extent in comets. Fireballs are due to multi-tons 
boulders some of which are believed to stem from comets as well. 

Spacecraft missions to comets must take careful precautions to survive the 
hazards of the cometary dust cloud. Special dust shields have been designed 
to protect the spacecraft. Nevertheless, the dust environment can have detri- 
mental effects on some aspects of the missions. The Giotto spacecraft flew by 
comet IP/Halley at a distance of 600 km and a speed of 70 km/s. During this 
close fly-by Giotto was hit by mm-sized dust particles. The impacts caused 
nutation of the spacecraft spin, and the data transmission was disrupted. 
Some experiment sensors suffered damage during this fly-by. Similarly, the 
350-kg impactor of the Deep Impact probe experienced attitude disturbances 
from dust grains before it hit the nucleus of comet 9P/Tempel 1. Owing to 
the low velocity of Rosetta relative to the comet, the consequences of dust 
impacts will be much less severe than for the fly-by missions. But detailed 
knowledge of the dust environment will be vital for the planning of space- 
craft operations at the comet and is thus of crucial importance to optimise 
Rosetta's scientific return. 

One of the problems in characterising the dust environment of a comet is 
that information on the nucleus, its dust and gas release is very limited. Be- 
fore 1986, observations of cometary dust was the domain of astronomers. The 
method of Finson and Probstein [T| was the first to be used to determine the 
size distribution of dust from observations of the tail at visible wavelengths. 
High resolution astronomical images of cometary comae revealed jets and 
other structures in the inner parts, some of which formed spirals which ro- 
tated like water from a lawn sprinkler indicating discrete dust emissions from 
localised active parts of the nucleus. A consequence of observing in visible 
light is that the results are biased to particle sizes in the range of 1 to 10 
/um, because much smaller and much larger particles do not contribute sig- 
nificantly to the scattered light [2]. With the extension of the observable 
spectral range to infrared wavelengths, also the thermal emission of dust be- 
came accessible to astronomers. It revealed information on the abundance 
of larger grains and on the mineralogical composition of the dust, the latter 
from characteristic spectral features in the near and mid-infrared range. 
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Breakthroughs in understanding cometary constituents came with the 
space missions to several comets: Giotto and two VeGa spacecraft to comet 
IP/Hahey in 1986, Deep Space 1 to comet 19P/Borelly in 1999, Stardust 
to comet 81P/Wild 2 in 2004, and Deep Impact to comet 9P/Tempel 1 in 
2005. Water and CO were identified as the main species in the gas, and dust 
particles made of carbonaceous and silicate materials ranging from nanome- 
tre to millimetre sizes were detected. Active areas and corresponding dust 
jets were identified in spatially resolved images of some of the comets visited 
by spacecraft [21 HI O For 67P/C-G, however, such detail will only be 
observed when in 2014 the Rosetta spacecraft reaches the comet. Until then 
any information on the dust environment has to be derived from astronomi- 
cal observations of the target comet or by assuming correspondence to other, 
better studied comets. 

The purpose of the present paper is twofold: on the one hand, to give an 
overview of the current knowledge of the dust environment of comet 67P / C- 
G, and on the other hand, to provide estimates of such quantities as the 
spatial density, flux, and speed of dust in the coma as functions of location 
and time. These values are meant to support the planning of measurements 
of the instruments on board Rosetta. Section [2] contains an overview of the 
available astronomical observations of comet 67P/C-^G, and measurements 
of the albedo and temperature of dust from 67P/C-G are presented. In addi- 
tion, the phase function, size distribution and radiation pressure efficiency of 
cometary dust are defined and discussed. In Section [3l several methods are 
presented to derive dust properties from modelling of images of the cometary 
tail or trail, and results for 67P/C-G that were obtained by several authors 
are discussed. Finally, in Section |4l the ESA Cometary Dust Environment 
Model [3 m m [To] is briefly described and results from applying it to comet 
67P/C-G are presented. 



2 Observations of 67P/C-G Dust and Dust Properties 

In this section, we give an overview of the available observational data con- 
taining information on the dust of 67P/C-G, and we introduce and discuss 
the major quantities that can be measured by means of such observations. 
Some of these quantities (the dust size distribution, the radiation pressure 
efficiency and the emission speeds) can - in the absence of in situ measure- 
ments - only be inferred through modelling of astronomical images, which is 
discussed in Section [3l Published astronomical data on the dust of 67P/C-G 
include images of the dust coma, tail and trail in both visible and infrared 
light. They are available from 1982 onwards. 
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2.1 Disambiguation: Dust tail, antitail, neckline, and 
trail 

In the following we give the definitions of some observational dust phenomena 
associated with comets as they are used in this paper: 

Outside the inner coma - in which the dust is accelerated by gas drag - 
dust dynamics is dominated by solar gravity and radiation pressure. Both 
forces follow a 1/r^-law (rh being the heliocentric distance) but act in op- 
posite directions. Consequently, their ratio /3 depends only on the material 
properties of the dust grains, such as size, composition, density, and shape. In 
general, radiation pressure is most efficient for particles of about the size of 
the dominant wavelength of the radiation. Such grains are driven away from 
the nucleus in the direction opposed to the Sun and trailing the nucleus, thus 
forming the comet's dust tail before they disperse into interplanetary space. 

The term antitail refers to a part of the tail that seems to point toward 
the Sun instead of away from it. Often, this is a projection effect that occurs 
when the observer is in such a position that part of the normal tail appears 
to be on the Sun-facing side of the nucleus. Viewed in three dimensions, there 
is no difference between an ordinary tail and a projection antitail, but the 
dust seen in the antitail tends to be the larger and older component of the 
dust in the tail. 

A neckline [HI [121 [E] consists of large particles emitted at a true anomaly 
of 180° before the observation: The orbital periods of large particles are 
similar to that of the parent comet. Their orbits are generally inclined with 
respect to the comet orbit, but the particles cross the orbital plane of the 
comet twice during each revolution around the Sun. One intersection point 
is the point of emission. The other lies on the line of nodes connecting the 
emission point and the Sun. The position of the second intersection point 
on the line of nodes depends on the emission velocity and (3 of the particle. 
Large particles emitted at a given time cross the orbital plane of the comet 
almost simultaneously, but at different positions along the nodal line. To an 
observer in - or close to - the comet orbital plane, they appear as a bright 
line, the neckline. Necklines can appear both in the Sun- and the anti-Sun 
direction i.e. can contribute to the tail as well as to the antitail. In the case 
of comet C1995 01 Hale-Bopp in early 1997 the neckline was mainly visible 
along the tail direction, but also gave rise to an antitail [Tl] . 

The dust trail of a comet consists of mm- to cm-sized particles that - 
because of low emission speeds and little sensitivity to radiation pressure 
- remain close to the comet orbit for many revolutions around the Sun and 
whose appearance reminds of an airplane contrail. Trails of eight short-period 
comets were first observed with IRAS in 1983 [TH[Tni[I7], one of them being 
that of 67P/C-G. 
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2.2 Morphology of Coma, Tail, Antitail and Trail 

The coma of 67P/C-G showed azimuthal brightness variations during both 
the 1996/97 and the 2002/03 apparitions [Ml [El |20l EJ. An example is 
displayed in Figured] The bright regions have been interpreted as border hnes 
of coma fans produced by active areas at different latitudes on the rotating 
nucleus. The 2003 fan pattern suggests the presence of 2 or 3 active regions on 
the nucleus. For the 2-activc-rcgion scenario, the rotation axis requires to be 
very much inclined to the orbital plane, while for the 3-active-region scenario 
a very wide range of rotation axis directions is possible (H. Bohnhardt, private 
communication) . The azimuthal direction of the bright features has been used 
to constrain the orientation of the rotation axis of the nucleus pTl [T5] . 




Fig. 1 Left: Broadband R image of comet 67P/C-G on 3 May 2003 obtained with 
FORSl at ESO/VLT/Ul [HjEOj. Right: Same image, but structurally enhanced by 
subtracting the azimuthally averaged coma profile from the measured brightness dis- 
tribution. Distinct features in the coma are clearly visible. The image is a reproduction 
of Figure 1 (lower left) in | 20 | , courtesy by R. Schulz. 



The tail of 67P/C-G was characterised by a thin, bright feature close to 
the projected comet orbit and pointing away from the nucleus in the direction 
opposed to the motion of the comet. This feature was first observed shortly 
after perihelion in August 2002 and prevailed at least until April 2006 when 
the comet had already passed aphelion [22l |23l EH EH [26] . Different interpre- 
tations of this phenomenon are discussed in Section|3]with the conclusion that 
it most probably was a very pronounced antitail due to the low inclination 
of the comet orbit with respect to the ecliptic. 

The dust trail of 67P/C-G was first observed with the Infrared Astro- 
nomical Satehite (IRAS) in 1983 [Til|Tll[nj with a reported length of 1.2° in 
mean anomaly and a width of 50000 km. In visible light, observations of the 
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67P/C-G trail were done in 2002/03 at heliocentric distances between 1.3 
AU and 3.1 AU ^ and in 2004 at rh = 4.7 AU (out-bound) [M]. Infrared 
observations of the 67P/C-G trail between 2004 and 2006 were obtained with 
the MIPS instrument of NASA's Spitzer Space Telescope at 24/Ltm [2711^ . 



(b) 



Fig. 2 Reproduction of Figure 1 in | 24 | . The images show the dust trail and necldine 
of 67P/C-G in April 2004. The data were obtained with the Wide Field Imager (WFI) 
at the ESO/MPG 2.2m telescope in La Silla with a total exposure time of 7.5 hours 
and without filter. The size of the images is 35'x4.7' each, corresponding to 1.1° 
in mean anomaly parallel to the orbit, (a) Unfiltered image, (b) Same image, each 
pixel being replaced by the average over a neighbourhood of 200 pixels (140") parallel 
and 10 pixels (7") perpendicular to the trail axis after removal of the nucleus. The 
filtering window is indicated in the upper right corner of (b). A detailed discussion 
of the data acquisition, processing and interpretation is given in |24| and modelling 
results are given in Section 13.31 



2.3 Albedo and Phase Function 



The geometric albedo p of an object is defined as the ratio of the intensities 
reflected backwards by the object and by a totally diffusely reflecting (i.e. 
Lambert scattering) disc of the same geometric cross section [5S], while the 
Bond albedo Ab is the ratio of the total light reflected from a sphere to 
total light incident on it The phase angle a is the angle between the 
directions of the observer and of the incident radiation as seen from the 
scattering particle. The phase function j(a) describes the ratio of intensity 
scattered in the a-direction to the intensity scattered at a = 0, and the 
integral of the phase function over 47r solid angle is called the phase integral, 
q. The Bond and geometric albedos are related by Ab = qp- 

To derive the dust phase function from astronomical observations, an ob- 
ject must be observed at different phase angles, i. e. at different times. Since 
the total cross section of dust is not constant with time, an appropriate nor- 
malisation is required, for which two methods are used. One method employs 
gas production rates measured simultaneously with Afp for normalisation 
[501 [5T1 132], assuming that the dust-to-gas ratio as well as the dust size 
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distribution and material properties remain constant over time. The other 
method, preferable but more laborious, normalises the scattered intensity to 
the simultaneously measured thermal infrared emission from the same vol- 
ume [31131133 [Ml ISZl [Si ISU- The general shape of the phase function of 
cometary dust is characterised by a distinct forward and a gentle backscat- 
tering peak and is rather flat at medium phase angles [40j . Divine |41j derived 
from data given in [37l [38l [39] the phase function shown in Figure [3| (solid 
line). The figure also shows the geometric phase function that describes the 
phase angle dependence of Afp and is discussed in Section [m 

Laboratory measurements and theoretical studies suggest that the dust 
albedo depends on particle size [42]. Earth based observations have so far 
not been suitable to investigate this dependence, because coma observations 
only provide data for the ensemble of particles of all sizes along a line of sight. 

The geometric albedo is derived from the simultaneous observation of the 
scattered visible and the thermally emitted infrared light, either directly at 
a = 0, or at multiple phase angles and assuming a given phase function. The 
geometric albedo of dust in the coma of 67P/C-G derived from the optical 
and infrared brightness was 0.04 at 1.25 ^m and 0.05 at 2.2 ^m [34], which is 
in accordance with a large sample of comets [43l [35] . From the low albedo it 
is inferred that there is no significant population of cold, bright (and possibly 
icy) grains that would contribute to the scattered light but not to the thermal 
emission |34| . There is, however, some indication that the geometric albedo is 
higher for comets beyond 3 AU |35| , which may point in the same direction. 
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Fig. 3 Solid line: phase function of an individual dust particle as given in [41| . but 
here normalised to j{a = 0) = 1. Dashed line: geometric phase function jgeo(a) [2] 
accounting for the anisotropy of an axis-symmetric coma with peak activity at the 
subsolar point. 
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2.4 Afrho and Gas Production 



The brightness of a cometary coma is proportional to the dust production 
rate. To infer the production rate from data obtained under different observa- 
tional circumstances, the measured brightness must be corrected for all other 
parameters on which it depends. The quantity Afp was defined for this pur- 
pose 131]. Afp stands for the product of albedo A = pj{a) (see Section [53]) , 
filling factor / of grains within the field of view, and the radius p of the 
aperture at the comet. It is measured as follows: 

Jsun ^ 

where rh and A are the heliocentric and geocentric distances of the comet 
during the observation, and p is the radius of the circular aperture on which 
the coma intensity /^us'T''' measured using a given filter, /iun is the in- 
tensity of the Sun at 1 AU heliocentric distance seen through the same filter. 
Provided that the dust particles move away from the nucleus on straight tra- 
jectories and are not subject to processes altering their scattering behaviour, 
Afp is independent of the employed aperture radius, of the heliocentric and 
geocentric distances, and - to the extent that the dust can be considered as 
"grey" - of the spectral band in which the observation was carried out. 

For an isotropic coma and discrete dust sizes Sj, Afp is related to the 
production rates Qd,j (sj) via the dust emission speeds Wd.j, the geometric 
albedo p, and the phase function j{a) [9]: 



A/p = 27rpj(«)E4— • (2) 



The relative magnitudes of the Qd,j rates are given by the size distribution. 

Afp depends on the phase angle of the observation, due both to the scat- 
tering properties of a single dust grain [32] and ~ unless the coma is isotropic 
- to projection effects [^. In a non- isotropic coma, the timescalc on which 
the particles leave a given field of view depends on the angle between the 
main emission direction and the line of sight, with higher measured Afp for a 
line of sight parallel or close to the main emission direction. The phase-angle 
dependence of Afp is then better described by a geometric phase function 
jgeoict) [S] than by that of a single particle. The geometric phase function is 
characteristic of the specific pattern of emission of a given nucleus. 

In practice, Afp often depends on the aperture size despite its definition 
[35] [TH] and references therein] . This implies that the brightness distribution 
in the coma deviates from the assumed l/p-profile. Possible causes for this 
deviation include changes in the physical properties of the grains as they 
travel outward (e. g. loss of volatiles or fragmentation), the action of radiation 
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pressure modifying the straight trajectories of smaU particles inside the field 
of view, or a long-lasting population of large particles [31] • 
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Fig. 4 Top left: observed Afp values as function of heliocentric distance rh- Top 
right: observed Afp values corrected for phase angle through division by the geometric 
phase function (Figure [3j. Power-law fit for phase angle corrected Afp: A/p(rh,a= 
0) = 2393 X (rh/1 AU)-5 °^ cm. Bottom: Measured H2O production rates with 
corresponding power-law fit QhjO (''h) = 3.4 x 10^® x (rh/1 AU)"^'^""^ molecules/s. 



Afp has been measured for 67P/C~G during three perihelion passages 
(1982/83, 1996/97 and 2002/03) gSlllillTlllilliinillllSni [IlE- Fig- 
ure [l](top left) shows all published Afp data as function of heliocentric dis- 
tance. The same values after correction for the phase angle dependence using 
the geometric phase function (Figure [3]) are displayed in Figure |4](top right) 
together with a power-law fit to the corrected values. Figure|4](bottom) shows 
observed production rates of H2O |34l EQI EU [52] and a power-law fit to the 
data. The exponents of the derived power laws are untypically steep com- 
pared with other comets, and more data of high quality will be needed to 
confirm them. Both Afp and the H2O production rate reach their maxima 



^ The data from these publications are summarised in the data base at 
http://berlinadmin.dlr.de/Missions/corot/caesp/comet_db.shtml. The site also in- 
cludes data obtained by amateur astronomers. 
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around 30 days after perihelion. No dust coma was detected beyond at least 
4.9 AU [Ml- 



2.5 Dust Temperature 

Assuming that a particle is characterised by the Bond albedo A-q at visible 
wavelengths and the emissivity e in the infrared, its temperature T at the 
heliocentric distance rh (in AU) is given by the equilibrium between absorbed 
solar and emitted thermal radiation: 

-2 =4e(Tr, (3) 

where a is the Stefan-Boltzmann constant and Iq = (1367 ± 2) W m~^ the 
solar flux at 1 AU |54j . The resulting equilibrium temperature is 

r(rh,AB,e) = 278.8 K fi^^V (4) 



A blackbody would be characterised by A-q = and e = 1. The temperatures 
of dust in the inner solar system are generally such that the main emission 
lies in the infrared. In practice, the temperature is derived from fitting a 
blackbody spectrum to measurements of the brightness at different infrared 
wavelengths. 

The coma of 67P/C-G was monitored at multiple wavelengths in the range 
of 1 to 20 /Ltm between September 1982 (1.50 AU preperihelion) and March 
1983 (1.88 AU postperihelion) [53]. The derived temperatures of the dust in 
the coma were throughout higher than those of a theoretical blackbody at the 
same heliocentric distance, which is generally attributed either to the presence 
of submicron-sized particles (i.e. smaller than the dominant wavelength range 
of the thermal emission) [SJ or to very porous aggregates of small grains 
|35[ 140] . For a given heliocentric distance, the colour temperature was higher 
post- than pre-perihelion, suggesting a change in the particle properties or 
in the dominant size. Additionally, an 8- 13-/im spectrum was taken on 23 
October 1983. It did not show a silicate feature, which is usually taken as an 
indication for the dominance of somewhat larger and more compact grains. 

Excess colour temperatures were also derived from IRAS observations of 
cometary trails at 12, 25, and 60 /im [T7]. For 67P/C-G, the derived temper- 
ature at 2.3 AU was approximately 14% above that of a blackbody at the 
same heliocentric distance, corresponding to e/(l — Ab) = 0.6 ± 0.2. 

Figure [5] shows the temperature measurements discussed above and - for 
comparison - the equilibrium temperature of a blackbody as a function of 
heliocentric distance. The data from the coma [33] and from the trail [T7] are 
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remarkably consistent given that coma and trail are generally assumed to be 
dominated by different particle populations (/im- versus mm-cm sized). 



Manner 



Blackbody 
al. (1985), preper. 
Manner et al. (1985), postper . 
Sykes and Walker (1992), postper. 
Sykes and Walker (1992), lower lira. 
Sykes and Walker (1992) , upper lira. 




Fig. 5 Colour temperatures of the dust from 67P/C-G derived from multi- wavelength 
infrared observations. The symbols "+" and "x" refer to observations of the coma 
before and after perihelion, respectively [34| . The asterisk indicates the temperature 
of the dust trail | 17 |. and the two dashed lines correspond - via Equation [4]- to the 
upper and lower limit of e/(l — As) derived from the same data. The solid line shows 
the temperature of a theoretical blackbody, described by Equation |4] with Ab = 
and e = 1. 



2.6 Dust Size Distribution 



The size distribution of eometary dust has been inferred from both astronom- 
ical images and in situ data. While the former yield a size distribution, the 
latter contain information on the masses of the particles. For optical images, 
the determined sizes scale directly with the particle albedo. The conversion 
from size to mass requires knowledge of the bulk density of the particles. 

A mass or size distribution can be specified in the form of either a differ- 
ential or a cumulative distribution. The cumulative mass distribution Fimo) 
gives information on the fraction N of particles that have a mass greater than 
some mass mo: 

N (m > Too) = F{mo). (5) 

The differential mass distribution /(to) characterises the relative abundance 
n of particles inside a mass interval [toi, m2]: 
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n{rni<m.<m.2) = j f {m) dm = F {rrii) — F {1712) ■ (6) 

mi 

If the mass of a particle can be converted to a size by a relation s(to), a 
corresponding differential size distribution g(s) exists: 

m2 s{m2) 

J /(m) dm = y" g (s) ds. (7) 

'"■1 s(mi) 

It is generally assumed that ~ at least in intervals of the total mass range 
covered by cometary dust - the distributions can be approximated by power 
laws. In the literature, both the exponents 7 of the cumulative mass distri- 
bution F{m) ^ m~'^ and a of the differential size distribution g(s) ^ s" are 
commonly used. For a constant bulk density in the concerned size interval we 
have m(s) oc s^ and dm cx s^ds. Hence the exponents 7 and a are related by 

a = -37-1. (8) 

The mass distribution at the nucleus is different from the one in the coma 
because of the size-dependence of the emission velocities. In general, large 
particles are more abundant in the coma than close to the surface because of 
their lower speeds. The relation between the size distributions in the coma 
and at the nucleus may be additionally complicated by fragmentation or 
evaporation of grains [SHI |311 [S5] and by an inhomogeneous distribution of 
surface activity [57] . 

In situ data on the dust mass distribution were obtained by the dust instru- 
ments on board the spacecraft VeGa 1 and 2 and Giotto at comet IP/Halley 
in 1986 (Sg EH SUEZ], on board Giotto at comet 26P/Grigg-Skjcllerup in 
1992 [SS], and on board Stardust at comet 81P/Wild 2 in 2004 [501 M\- 
The measured quantity is not the mass distribution of dust as released from 
the nucleus but the cumulative flux or fluence on the concerned instrument. 
The fluence is defined as the flux integrated over the spacecraft trajectory, 
it represents therefore an average mass distribution. The flux (or the fluence 
measured in only sections of the trajectory) showed signiflcant variation with 
time during both the IP/Halley and the 81P/Wild 2 fly-bys. 

Table [T] lists mass- or size-distribution exponents given in the literature. 
Still no general agreement has been reached on the interpretation of the data 
with respect to the dust mass distribution at the nucleus, but the authors 
listed in Table [1] do agree that the mass of dust in the coma is dominated 
by the largest emitted particles. The cumulative fluences registered at the 
various spacecraft show different exponents for large and small particles. The 
interpretations of this observation are not unanimous, either. 
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Table 1 Exponents of the cumulative fluence on a spacecraft, 7, the cumulative mass 
distribution, 7, and the differential size distribution, a, compiled from the literature. 
Bold values are taken directly from the publications, the remaining values were de- 
rived assuming the following relations. To translate 7 to 7, it is assumed that the 
speeds of small particles are limited by the speed of the gas and therefore indepen- 
dent of size. Hence 7 = 7 for small particles (indicated by "f")- For larger particles, 
the relation v oc m~^^^ is used, which was derived from a radially symmetric coma 
model |43| . resulting in 7 = 7 -|- 1/6 for large particles (indicated by "*"). The rela- 
tion between the differential size distribution and the cumulative mass distribution 
is given by a = —37 — 1. Results from modelling of astronomical images are not 
included here. For 67P/C-G a compilation of size distributions derived by various 
authors is shown in Figure [7](bottom). 
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10-19. 


.. 2x10-14 kg 





18 


0.18t 


-1.54 


m 






2x10-" . 


.. 10-" kg 





94 


1.11* 


-4.32 






VeGa 1 


10-19. 


.. 10-1=^ kg 





26 


0.26t 


-1.78 








10-12. 


. . 10-9 kg 


1 


19 


1.36* 


-5.07 






VeGa 2 


10-19. 


. . 1.6x10-13 kg 





26 


0.26t 


-1.78 








1.6x10-13 . 


. . 10-9 kg 





90 


1.07* 


-4.2 






DIDSY 


> 


20 lira 








-3.5 ± 0.2 


m 


26P/G-S 


DIDSY/GRE 


10-9 . 


.. lO-'' kg 





Q-+0.13 
^ ' -0.20 


0.44* 


-2.3 


m 


81P/Wild 2 DPMI 


10-14. 


.. 10-9 kg 





85 ± 0.05 


1.02* 


-4.05 ± 0.15 


[56] 




Stardust Samples 


10-1^ 


.. 10-3 kg 





57 


0.74* 


-3.21 


m 



2. 7 Radiation Pressure 



Both solar gravity and the radiation pressure force are inversely proportional 
to the heliocentric distance squared, and point radially away from the Sun. 
Radiation pressure can, therefore, be included in the equation of motion of 
a particle by introducing a modified potential substituting the gravitational 
constant, G, by G = G(l — /3), where /3 stands for the ratio of solar gravity 
to the radiation pressure force. It depends only on material properties of the 
dust grains, not on their distance from the Sun: 

B = 9pL ln\ 

' 16 ttcGMq ps' ^ ' 

Lq and Mq are the luminosity and mass of the Sun, and c is the speed of 
light. The grain has the bulk density p and the effective radius s, and it is 
characterised by the radiation pressure efficiency Qpr [Ml US] . Qpr represents 
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the absorption and scattering properties of the grain averaged over the solar 
spectrum. For homogeneous spheres and some other simple shapes, Qpy can 
be calculated in an exact way from the complex refractive index of the ma- 
terial [531 [53]. For more complicated structures, a variety of theoretical and 
experimental approaches exist to obtain Qpr [65]. While for particle sizes on 
the order of 0.01 to ly^tm, f3 depends sensitively on material, shape, struc- 
ture, surface properties, and size of the particles, it is approximately constant 
for much smaller particles, and proportional to l/(ps) for large ones, i.e. in 
the geometric optics regime. This implies that for sub-millimetre and larger 
particles, Qpr is to first order independent of size [62l |64] . 



3 Results of Modelling the 67P Dust Tail and Trail 

In this Section, we give an overview of the efforts undertaken to derive prop- 
erties of the 67P/C-G dust from modelling astronomical observations of the 
thin, bright feature described in Section [221 variantly referred to as neckline, 
antitail, or trail. The main quantities to be constrained by the models are 
the dust size distribution, the radiation pressure efficiency, and the emission 
speeds. The derived values will serve as input to the hydrodynamic coma 
model used in Section [4| to predict the dust environment inside the coma. In 
the following, we summarise the results from neckline photometry applied to 
images taken in 2002 /03 [22l [23| , from analysis of the antitail as observed in 
May 2003 [5^ , and from simulations of the dust trail in visible and infrared 
wide-field observations between 2002 and 2006 [25 l [26 l [27] . 

The position of a dust particle in a cometary tail is a function of the radi- 
ation pressure coefficient /?, the emission time, and the emission speed; and f3 
is closely related to the mass of a particle. Hence, the effect of solar radiation 
pressure is similar to that of a mass spectrometer. To understand the forma- 
tion of dust tails, the concept of synchrones and syndynes was introduced by 
Bredikhin [67] . Synchrones are the positions of particles of different (3 emitted 
at a given time, while syndynes describe the positions of particles of fixed /3 
and varying emission time. Both terms refer to hypothetical particles emitted 
with zero velocity relative to the nucleus. Since, realistically, the initial veloc- 
ity of a dust grain is different from zero, the resulting synchrone or syndyne 
will have a finite cross section with a radius proportional to the product of 
dust emission speed and time elapsed since emission. In particular, grains 
released at a given time with isotropic speed will form a spherical shell the 
centre of which moves along the appropriate syndyne. 

This description was used by Finson and Probstein [T] [55] to derive prop- 
erties of the dust size distribution from the brightness patterns observed in 
cometary dust tails. Their method is limited to recently emitted and small 
particles, because it neglects tidal effects and others (e.g. direction-dependent 
emission speed or production rate) that cause a dust shell to divert from the 
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spherical shape [TTl IBH [13] . Various approaches to surmount these hmitations 
are described in the fohowing. 



3.1 Neckline Photometry 

Fulle et al. [22] analyse 67P/C-G images obtained in March 2003 (rh = 2.6 
AU) and before, when the comet was active and exhibited a significant coma. 
They employ both an analytical theory of neckline photometry and an inverse 
Monte Carlo model. 

The former is an analytical method to infer the emission speeds and abun- 
dance of particles as a function of the radiation pressure parameter /3 |70| . It 
is applied to an image of the comet obtained with the Schmidt Telescope at 
the Thiiringer Landessternwarte (TLS) in Tautenburg on 27-28 March 2003. 
The method relies on the assumption that the bright narrow feature in the 
tail is a neckline, and that all intensity observed in this feature is due to 
dust emitted within an interval of 10 hours at a true anomaly of 180° before 
that of the observation, which corresponds to 5 May 2002 (105 days before 
perihelion) and rh = 1.8 AU. The authors find that - in this particular time 
interval - both the mass and the cross section of dust emitted by the comet 
were dominated by particles in the size range of 1.5 to 10 mm, and that 
the exponent of the differential size distribution, a, was between -3.5 and 
-3. From the fact that the described feature had been constantly observable 
since perihelion in August 2002, they infer that the characteristics of the dust 
emission as recorded in the neckline did not change between 3.6 and 1.7 AU 
before perihelion. 

The second method - the inverse Monte Carlo model - consists in a least- 
squares fit of simulated images to the measured surface intensity |71j . This 
method yields a set of time-dependent dust parameters (size distribution, 
emission velocity, size range, Afp, and dust production rate). The obtained 
solution is unique in the mathematical sense of a least-squares fit; its physical 
probability remains harder to evaluate. The method is applied to an image 
obtained with the Telescopio Nazionale Galileo (TNG) on La Palma during 
the same night as the TLS image. Quantitative interpretation of the model is 
achieved through matching of the modelled Afp to the value measured during 
a different observation close in time [49]. The authors [22] conclude that 
the time dependence of all parameters must be asymmetric with respect to 
perihelion in order to match the observation. The size distribution exponent 
drops from —3.5 before to —4.5 after perihelion, the mass loss rate from 
200kg/s to lOkg/s. The emission speeds decrease by a factor of three during 
this time interval (from 3m/s to Im/s for a particle of 1mm radius), and 
the dominant size by two orders of magnitude. Fulle et al. [22] suggest the 
existence of two different source regions: one emitting large particles and only 
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active before perihelion, and one emitting mostly small particles becoming 
active after perihelion. 

A decreasing production rate was also found by Moreno et al. [23] applying 
a similar method to the same image after calibration. 



3.2 Antitail Analysis in the Hypersonic Approximation 

The dust size distribution has been evaluated by analysis of the brightness 
profile along the 67P/C-G tail using the Finson-Probstein model in the hyper- 
sonic approximation |66| , which was applied to an image (Figure [1]) obtained 
with ESO/VLT on 3 May 2003 

At the time of the observation, 67P/C-G was at a heliocentric distance of 
2.9 AU postperihelion. A plot of synchrones and syndynes (Figure [S]) shows 
that all synchrones corresponding to ejection before 1.5 AU postperihelion 
(23 October 2002) appear projected in the direction towards the Sun, forming 
an antitail. They all are - in projection - inclined against the neckline by less 
than 0.8°, which places doubts upon the assumption that all dust in the 
bright feature was emitted during an interval of merely 10 hours, on which 
the first model described in Section [Ol relies j22j . 

The hypersonic approximation holds when the extent of the dust shells 
along a synchrone is small, such that a unique relation between the position 



[lOOOkm] [1000km] 




Fig. 6 For 3 May 2003 synchrones are shown projected in the comet orbit plane 
(top) and image plane (bottom). The solid lines correspond to synchrones ejected 
in steps of 10° in true anomaly starting 5° in true anomaly before the observation 
time back to an ejection time that corresponds to a heliocentric distance of 3 AU 
before perihelion. Synchrones corresponding to ejection times after mid January 2003 
appear in the direction opposed to the Sun. The synchrones of particles ejected before 
mid December 2003 appear in Sun direction. Only the synchrones of particles ejected 
in between (shown as long dashed) have significant angles with the Sun or anti-sun 
direction. Projected in the comet orbit plane (top panel) also syndynes are shown 
(dashed). The /3-values corresponding to the syndynes are annotated for the larger 
particles. The /3 values of the smaller particles are LSxlO"**, 3x10"'', and 1x10"'^. 
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on the synchronc and /? can be assumed Q]. It is then possible to infer the 
size distribution of the particles ejected from the nucleus from the brightness 
profile along the synchrone, if one of the following two conditions is met: 
(1) Either the antitail originates from the superposition of many synchrones 
that - in projection - have only a small offset from each other compared 
with the extension of the dust shells perpendicular to the synchrone. Or (2) 
the brightness along one direction in the image plane is dominated by one 
synchrone alone. A generalisation of the original formulation of the hyper- 
sonic approximation [T] - which applied only to spherical shells - has been 
developed for tidally distorted elliptical shells, assuming radially symmetric 
particle emission [66] . 

The intensity along the antitail on 3 May 2003 could be fitted by a power 
law with an exponent of —0.4, which was translated to a = —4.1 for the 
differential size distribution between at least 3 AU before and 1.5 AU af- 
ter perihelion in 2002 |66j . To put this result in a wider context, we have 
compared the measured brightness distribution along the antitails of differ- 
ent comets rather than the derived size distributions, because we found that 
different authors use different models to derive the size distribution from the 
measured brightness exponent. We found that the brightness variation mea- 
sured for the antitail of 67P/C-G is rather typical compared with the antitails 
of other comets, such that also the dust size distribution of 67P/C-G may 
be typical. 



3.3 Trail Analyses 

Several wide field images of the dust along the orbit of 67P/C-G were ob- 
tained between 2002 and 2006 in both optical and infrared light, and their 
analyses are summarised in the following. 

Three images - one optical and two infrared at 24 ^m - with fields of view 
of about half a degree were analysed using a generalised Finson-Probstein 
approach that takes into account the tidal deformation of dust shells on long 
time-scales [25] . The optical image was obtained in April 2004 with the Wide 
Field Imager at the ESO/MPG 2.2m telescope on La Silla at rh = 4.7 AU 
(Figure [5]) , and the infrared images were taken by the MIPS instrument on 
board the Spitzer Space Telescope of NASA in August 2005 and April 2006 
(rh= 5.7 AU in both cases). All three images were thus taken when the 
comet was close to aphelion and not active, such that no particles larger 
than approximately 100 /im are expected to be present in the fields of view. 
Simulated images were generated taking into account dust emitted during all 
seven perihelion passages since the last close encounter with Jupiter in 1959. 
The time-dependence of the dust production was modelled on the observed 
time-dependence of Afp. The emission speeds and surface activity were as- 
sumed to be isotropic, and the relative dependence of the emission speeds 
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on size and heliocentric distance was obtained by help of the hydrodynamic 
coma model described in Section 0] of this article, and from the observed 
time-dependence of the water production rate. Within this framework, the 
parameter values found most suitable to reproduce both the brightness pro- 
file along the trail and its width were the following: The emission speeds of 
a 1mm radius particle at perihelion ranged from 6m/s to 12m/s, decreasing 
by a factor of 10 until 3 AU. The size dependence of (3 could be characterised 
by Equation[9]with Q-pxl p between 1 and 3 cm^/g, which - assuming Qpr = 1 
for large particles - corresponds to bulk densities between 0.3 and 0.9 g/cm"^. 
The size distribution exponent for particles larger than 100 /xm was in the 
range of -4.1 to -3.9. In order to reproduce both the optical and the infrared 
images, a geometric albedo for visible light of p = 4% was required if the 
dust was assumed to be emitting as a blackbody (e = 1 in Equation H} , but 
•p = 10% if e/ (1 — A) ~ 0.6 [331 [T7]. Lower limits for the production rates of 
particles with s > 100 /im ranged from 100 kg/s at perihelion to 0.2 kg/s at 
3 AU, corresponding to Afp values of 4 and 0.05 m, respectively. This implies 
that particles larger than 100 /im would have contributed at least 50% of the 
total Afp observed while the comet was in the inner solar system, which is 
difficult to reconcile with the size distribution exponent on the order of -4. 

Three optical images (R-band) were obtained between September 2002 and 
February 2003 with the 1.05-m Schmidt telescope of the Kiso Observatory 
at Nagano, Japan [26]. The difference in position angle of freshly emitted 
dust and the trail (dust from previous apparitions) was 1° or larger in these 
images, thus - in contrast to images taken at later dates - a clear separation 
of both dust populations is possible. Taking into account dust emitted after 
aphelion in 1986 and assuming p = 4%, p = 1 g/cm'^, and Qpr = 1, the 
images were best reproduced by a model with cone-shaped emission with a 
half opening angle between 45° and 90°, a size distribution exponent of -3.5 
with particles in the range between 6 pm and 5 mm, and dust production 
rates of about 200 kg/s at perihelion and 15 kg/s at 3 AU. The emission 
speeds were assumed to be proportional to s~^/^ and absolute values 

for a 1 mm particle ranging from 8 to 18 m/s at perihelion and from 5 to 12 
m/s at 3 AU. 

Two further observations were analysed by help of a Monte Carlo model 
[27] : an optical image (Gunn r' filter) obtained in June 2003 with the 5-m 
Hale telescope at Palomar Observatory, and an infrared (24 /im) image made 
with Spitzer/MIPS in February 2004. The trail of particles from previous 
perihelion passages was not detected in the Palomar image, but was visible 
to Spitzer. The applied model includes the assumptions of Qpr — 1, p ~ 1 
g/cm^, dust production rates Qd oc r^^'^, and emission speeds v oc \/pJr\. 
Particles larger than 0.5 pm and emitted after March 1993 were included 
in the analysis. The images were best fitted by a cosine-shaped distribution 
of the surface activity, peaking at the subsolar point, and a differential size 
distribution exponent of -3.5. The emission speeds of a 1mm particle at the 
subsolar point varied between 10 m/s at perihelion and 7 m/s at 3 AU. 
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3.4 Summary of Modelling Results 



Figure [7] summarises the dust production rates, emission speeds and size 
distributions derived by help of the various models described in Sections 13.11 
to l3.3l We wish to stress that - with the exception of the TNG image obtained 
in March 2003 - each model has been applied to a different set of images, 
and each model was able to reproduce the images that were analysed by it. 
However, the ranges of parameter values derived from the different models are 
considerable, such that it is at present difficult to derive a consistent picture 
of the CG dust environment on the basis of these results. Future work on this 
matter should focus on finding a set of parameters that is able to reproduce 
all available observations of the CG tail and trail. 




-3OD -20a -100 100 200 3DD -3DD -200 -100 100 200 300 




-300 -2O0 -100 100 200 300 



Fig. 7 Comparison of dust characteristics derived from the different models described 
in Section 121 Two lines of the same style represent upper and lower limits. Top left: 
dust production rates. Top right: emission speeds of a dust grain with 1 mm radius. 
For anisotropic models, the peak emission speed (assumed above the subsolar point) 
is given. Bottom: Exponents of the differential size distribution. 
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4 Dust Environment Model and Predictions for 
67P/C-G 

In this section, the state of dust in the coma of 67P/C-G is predicted by 
means of a hydrodynamic model [3 [8j [9l [10] . Three emission scenarios and 
three size distribution exponents are studied, and estimates of the fluxes on 
a spacecraft for several types of trajectory and the radiation environment in 
the coma are presented. 



4-1 Coma Model and Parameters 

The employed model can be considered as the simplest possible physically 
consistent model of the inner coma: The comet nucleus is assumed to be of 
spherical shape, and the activity distribution on its surface is axis-symmetric 
with respect to the comet-Sun line. The flows of gas and dust around the 
comet are calculated rendering at any position in the coma the gas density, 
velocity, and temperature, and the number density and velocity of dust par- 
ticles of different mass classes. Each class represents particles having masses 
within one decade in the range of 10-20 to lO^'kg. All dust particles of one 
class are assigned the same representative mass and radius. It is assumed that 
the particles have spherical shapes and a constant bulk density of 1000 kg/m"^. 

The numerical method used to compute the dynamics of gas and dust 
takes advantage of the fact that the influence of the dust on the gas flow 
is of minor importance. In a first step, the gas flow is calculated without 
taking into account the presence of dust. In a second step, dust trajectories 
are integrated considering the gravitation of the nucleus and the gas drag 
force. 

The gas flow is calculated under the assumptions that the gas is in thermal 
equilibrium everywhere, and that the mean molecular mass is constant across 
the coma. It is given by the mean of the masses of the two most common 
molecules, II2O and CO, weighted by their overall abundances in the coma. 
The production rates of these species are input parameters to the model. 
To first order, the CO-activity can be treated as independent of heliocentric 
distance, and the constant value of 10^^ molecules/s is adopted for 67P/C-G. 
An upper limit for the CO-production of 10^^ molecules/s at 3 AU has been 
inferred from radio observations [75]. The water activity is assumed to vary 
with heliocentric distance and time as indicated by observations (Figure [U 
bottom). The dust activity is scaled such that the observed Afp is reproduced 
for the employed size distribution and albedo. 

Since, at present, the size distribution estimates vary considerably (Sec- 
tion [33]), we consider in the following the same range of mass distributions as 
in earlier models for 46P/Wirtanen [7]: two "extreme" and a "nominal" one. 
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All are characterised by different exponents for light and heavy particles, and 
the analytical form used for the cumulative mass distribution is as follows 

This function is specified by the positive parameters a, b, c, and mt- The 
exponent —7 of the cumulative mass distribution tends towards —a for heavy 
particles (to <^ rrit) and towards ~ab for light ones, mt is the mass where 
the transition between the two exponents takes place, and c determines the 
sharpness of the transition. The three mass distributions considered in this 
paper are different only by their slopes for heavy particles. For light particles 
we use 7 = afe = 0.26 throughout (corresponding to a = —1.8), for the 
transition mass mt = 10^^'^ kg, and c = 2. The slopes for large particles are 
as follows: the (velocity-corrected) fit to the fluence measured on Giotto [58] 
gave a ~ —4.3 (7 ~ 1.1) which we use as one of the "extreme" distributions. 
From the same data, an exponent in the range of a G [—3.7, 3.3] has been 
derived [57|, wherefore we adopt a ~ —3.7 (7 = 0.9) as the "nominal" and 
a = —3.3 (7 = 0.8) as the other "extreme" distribution. It must, however, 
be emphasised that for a > —3.5, the observed brightness is dominated by 
light scattered by large particles [7] . Figure [5] show different moments of the 
mass distribution. 

In the following, model predictions are given for several heliocentric dis- 
tances: the perihelion distances in 2015 (rh = 1.24 AU) and 2002 (rh = 
1.29 AU), and 1.87 AU, 2.5 AU, and 3 AU. The values used for the water pro- 
duction rate, for Afp, and for the phase angle a are listed in Table [2l We 
use a geometric albedo of 4% and the phase function shown in Figure |31 
While the production of CO is assumed to be isotropic, for the release of 
H2O three different surface activity distributions are taken into account. The 
first represents a homogeneously composed surface such that the strength of 
gas production is determined by the amount of solar energy received by a 
given surface element (insolation driven model). In the second scenario, a jet 
is simulated describing the surface activity by a Gaussian profile centred at 
the subsolar point [57l [73] . This model is expected to yield upper limits for 
local quantities inside a jet, while the insolation driven model gives estimates 



Table 2 Water production rates Qh,Oi A/p-values and phase angles a for those 
heliocentric distances rh used to obtain the model predictions in this section. 



r-h [AU] 


Qh,o [molecules/s] 


Afp (a) 


[cm] a[c 


1.24 


1.3 X 10^8 


400 


35 


1.29 


1.0 X 10^** 


300 


35 


1.87 


1.2 X lO^'^ 


65 


32 


2.5 


1.8 X 10^'' 


20 


5 


3.0 


6.4 X 10^5 


5.8 
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Fig. 8 Different moments of the mass distributions studied in Section |4l Top left: 
number fraction. Top right: area fraction. Bottom: mass fraction. To convert mass to 
radius, a bulk density of p = 1000 kg/m'^ was assumed. 



for quantities averaged over time and space. In addition, an even simpler 
radially symmetric model is used in which the water activity is distributed 
equally over the comet surface. 

The average radius of the nucleus of 67P/C-G as given in the literature 
is in the range of 1.7. .3. 2 km [H [751 ESI EZl El EH [27]. The geometric 
albedo of the nucleus in R-band is in the range of 0.045 to 0.06, and the bulk 
density of the nucleus is estimated to 370kg/m^ [T^- The following results 
were obtained for a spherical nucleus having a radius of 2 km and a geometric 
albedo of 4%. 



4-2 Dust Terminal Speeds 

The terminal speeds of dust as a function of particle size for the three different 
emission scenarios are shown in Figure [9] For comparison, the corresponding 
values for comet 46P/Wirtanen, the former Rosetta target, are shown as well. 
Note that all curves are at perihelion, hence the comets are compared at 
different heliocentric distances. While the terminal speeds for 46P/Wirtanen 
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Particle radius [m] for p = 1000 kg/m 
10"' 10"* 10"^ 10"' 10"^ 10"^ 10" 



Wirtanen, r.s. 

CG, r.s. 
CG, i.d. 

CG, jet 



0.1 ^ ^ ^ 

10 ^ 10 * 10 10 ^ 10 ^ 10° 10^ 10^ 10'' 

Mass-to-cross-section ratio [kg/m^] 

Fig. 9 Terminal speeds as functions of mass-to-cross-section ratio for various emis- 
sion scenarios (r.s.: radially symmetric, i.d.: insolation driven). The corresponding 
particle radii for a bulk density of p = 1000 kg/m'' are indicated at the upper margin. 
Values for 67P/C-G are given at 20 km from the nucleus centre above the subsolar 
point and at perihelion (rh = 1.29 AU). For 46P/ Wirtanen they are at 10 km from 
the nucleus centre and at the perihelion distance of rh = 1.06 AU. 



are given at 10 km from the nucleus centre, those for 67P/C-G are at 20 km, 
because the dust particles have not quite reached their terminal speeds at 
lower heights. 

In the radially symmetric models, the terminal speeds are higher in the 
case of 46P/Wirtanen than of 67P/C-G. The reason is that both comets have 
similar gas production rates although 46P/Wirtanen is smaller. Hence the gas 
directly over the surface is denser and the dust reaches higher velocities. The 
ratio between the speeds expected for 67P/C-G and 46P/Wirtanen in the 
radially symmetric model ranges from 0.8 for small particles to 0.5 for large 
ones. The speeds shown in Figure IH] were computed for a dust bulk density of 
p = 1000 kg/m'^. Scaling to other particle densities can to first order be done 
assuming that particles with the same cross section to mass ratio reach the 
same velocities. For the larger particles, a power law of the following form 
was fitted to the calculated speeds: 

^5)=«oQ (11) 

with So = Imni. The resulting parameters for the different scenarios are 
given in Table [31 For the fits, all dust classes with particle radii s > 0.1 mm 
have been used. The class with the largest liftable particles was not in- 
cluded because of the strong influence of the gravitation of the comet. For 
46P/Wirtanen, the exponent d is closer to the expected value of 0.5 than for 
67P/C-G because of the latter's larger nucleus and resulting stronger gravity. 
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Table 3 Results of fitting a power law to the dust terminal speeds for large particles 
above the subsolar point at perihelion. Values are given for different emission models: 
radially symmetric (r.s.), insolation driven (i.d.), and jet model. 



Model rh [AU] vq [m/s] Exponent d 



Wirtanen r.s. 


f.06 


7.0 


0.51 


C-G r.s. 


f.29 


3.9 


0.53 


C-G i.d. 


1.29 


8.0 


0.53 


C-G jet 


1.29 


22.2 


0.52 



4.3 Dust Flux at 1.29 AU 



Here we give numerical values for the fluxes of gas and dust on a surface 
directed towards the nucleus at 20 km from the nucleus centre above the 
subsolar point for various mass distributions and emission scenarios. As the 
speeds of gas and dust outside the acceleration zone are almost constant, the 
fluxes there can be scaled by an inverse square law. Note that the values given 
for the gas are only valid as long as the gas is cold, i.e. far from the nucleus. 
Within a few comet radii from the surface, the velocity distribution of the 
gas should be taken into account. Consequently, there will also be a flux on 
surfaces which arc directed away from the flux direction. The numerical values 
are given in the appendix of }66j and visualised in Figure 1101 The number 
and mass fluxes (upper panels Figure llOp in are relevant for different sensors 
of the GIADA instrument. The lower left panel shows the surface coverage 
by dust which is important for estimating thermal effects on spacecraft and 
instruments. The lower right panel describes the momentum transfer onto 
the spacecraft by the dust. 



4.4 Dust Collection along Sample Spacecraft 
Trajectories 

During most of the mission, regions of strong dust emission in the coma 
will be avoided, because impacts of dust particles constitute a hazard to 
the Rosetta spacecraft. This hazard includes damage of sensitive surfaces 
by cratering and erosion due to high velocity dust impacts, contamination 
of surfaces and modification of their optical and thermal properties due to 
sticking of dust, and even mechanical blockage of gears and hinges by larger 
grains. This hazard has been taken into account in the design of the Rosetta 
spacecraft and its instruments, e.g. by painting the spacecraft black so that 
sticking dust will not significantly change the thermal properties of the space- 
craft, and by employing movable shutters in front of sensitive surfaces such 
as cameras. Nevertheless, in the vicinity of the comet, spacecraft operations 
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Fig. 10 Fluxes over particle mass for various mass distributions and emission ge- 
ometries. Top left: number flux, top right: mass flux, bottom left: area flux, bottom 
right: momentum flux. 

will keep exposure to dust to a minimum [79]. In order to fulfil the dust col- 
lection requirements of the COSIMA and MIDAS instruments the spacecraft 
occasionally will have to pass through regions of high dust density. The de- 
sign of these passages will depend on the collection requirements of the dust 
instruments. 

For simple dust collection trajectories, the surface coverage by dust (as- 
suming all dust sticks to the surface) has been calculated for the different 
dust emission scenarios and mass distributions described in Section [4. II The 
considered trajectory types are one revolution on a circular orbit, a parabolic 
orbit, and a straight line trajectory passing the nucleus at a speed v = 1 
m/s. All trajectories are assumed to pass at closest approach (rcA = 4 km 
from the nucleus centre) over the subsolar point, i.e. the region of highest 
dust density in our model. The fraction of the surface covered by dust after 
one passage of the nucleus is given in Table H) For other closest approach dis- 
tances the surface coverage levels scale with for Keplerian orbits (circles 
and parabola) and with and for straight line passages. The values 
in Table |4] refer to different dust production rates, i.e. different heliocentric 
distances. 

A surface coverage of one implies roughly a mono-layer of differently sized 
dust grains, most of which are grains in the 1 to 10 /zm size range (cf. FigurelU 
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top right). Since the dust instruments arc sensitive to different sizes of dust 
grains it is useful to determine the dust coverage by particles above a certain 
minimum size. This is achieved by applying the scaling factors for different 
minimum dust sizes given in Table [S] Useful surface coverage levels of dust 
collectors range from 10"'* to 10-2, i.e. a 1 cm collection surface will contain 
approximately 100 to IC* particles of 10 /im in size. According to Table[5]these 
particles constitute about 50% of the covered surface (for the size distribution 
with a ~ —3.7). The other half of the covered surface consists of more finely 
dispersed smaller grains. 

The COSIMA instrument is sensitive to dust grains of 10 ^m and has 
23 individual collectors. The MIDAS instrument is sensitive to submicron- 
sized dust grains and has 60 individual collectors. The exposure of these 
collectors to the cometary dust flux and the collection of sufficient numbers 
of dust particles will at the same time lead to significant dust coverage of 
all spacecraft surfaces facing the comet. During one dust collection passage 
of Rosctta at large heliocentric distances (2.5 and 3 AU) an individual dust 
collector may reach only a total dust coverage of approximately 2 x lO"'^ and 

Table 4 Total surface coverage of the nucleus-directed side of the spacecraft due to 
dust particles of all sizes. The table lists the surface coverage for different heliocentric 
distances and different kinds of orbits. At heliocentric distances of 2.5 and 3 AU 
the water production rate is comparable with the background CO production rate, 
therefore, dust emission in form of a jet is not considered. We expect that at these 
heliocentric distances particles on the order of 100 ^tm can still leave the comet. All 
orbits have the point of closest approach to the nucleus over the subsolar point with 
a distance rcA = 4 km from the nucleus centre. The considered types of orbits are a 
parabola, a straight line at constant speed of 1 m/s, and a circle. For the parabola 
and the circle only the part of the trajectory over the comet day side is taken into 
account. 
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Table 5 Fraction of area flux due to particles greater than mass mmin or radius s 
for the different mass distributions. 
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7 X IQ-'*, respectively. Close to perihelion, however, the comet-facing side of 
the spacecraft will be almost completely covered by dust during a single dust 
collection passage. 



4-5 Radiation Environment 

In the following we predict the radiation levels received in the vicinity of 
the comet nucleus due to the presence of dust. It is beyond the scope of the 
present work to introduce a detailed wavelength-dependent model. Instead, 
values will only be given for the total amount of scattered visible light and 
thermally emitted radiation. The crucial parameters for these estimates are 
the albedo, phase function and temperature of the dust particles. We adopt 
the phase function derived for comet IP/Halley [41] as shown in Figure [3l a 
dust geometric albedo of 4% [34], and a temperature at perihelion of 285 K 
(see Figure [S]). 

A priori it is not clear that it is valid to assume that the temperature of a 
particle of a given size is independent of its position in the coma, because the 
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gas molecules impacting on the particle surface give rise to a heat exchange 
between the gas and the dust phase. This effect was taken into account in 
early models [SHI IHI] , but was found to be of minor importance later [331 . 
Neglecting the heat exchange, the temperature of a dust particle is given 
by the equilibrium between incident radiation and thermal emission of the 
particle. If the coma is optically thin, which is shown to be true for the dust 
coma of 67P/C-G below, the incident radiation is dominated by solar illumi- 
nation. Hence it is justified to treat the particle temperature as independent 
of the position in the coma. Assuming furthermore that the physical prop- 
erties of the particles on their way from the inner to the outer coma do not 
change significantly, the infrared spectrum of dust near the nucleus can be 
approximated by the spectrum observed from Earth. 

The thermal emission spectrum of cometary dust particles can only ap- 
proximately be represented by the spectrum of a blackbody, because particles 
emit efficiently only at wavelengths smaller than their size. Consequently, the 
smallest (or very porous) particles have temperatures much higher than a 
blackbody at the given heliocentric distance, and it is not physically consis- 
tent to set the dust temperature to a constant value independent of particle 
size. However, for the purpose of estimating the amount of radiation received 
by an observer, only the properties of the dust particles as an ensemble are 
of importance. The relative contribution to the total cross section of particles 
of a given size in most of the coma volume is the same as in the ensemble 
seen by an Earth-based observer. This was found with the present model 
which consistently models the dynamics of particles of different sizes. To cal- 
culate the intensity received inside the coma it is therefore valid to assign 
every particle independently of its size the ensemble properties as measured 
in Earth-based observations, which is henceforth done for the optical and 
thermal characteristics of the particles. 

Compared with the results given in the previous subsections - which were 
strongly controlled by model parameters only very indirectly accessible to 
ground-based observations (e.g. particle size distribution, particle density) 
- the results presented here only depend on parameters readily measured 
from ground. For example the estimates of the scattered visual radiation are 
governed solely by the Afp value and the dust phase function. It may, at first 
glance, be surprising that the results are not influenced by the dust albedo. 
However, because in the present model the overall dust activity is adjusted in 
order to match an observed Afp value, the predictions of the scattered visual 
radiation inside the coma are merely a scaling of the Earth-based observations 
to a different geometry. To give an idea of the accuracy of the results obtained 
for the scattered radiation: The present model applied to comet IP/Halley 
during the fly-by of Giotto reproduces the measurements of the Halley Optical 
Probe experiment [83l|84] by a factor of only 1.5. By contrast, the predictions 
of the optical thickness and of the thermal radiation by the dust depend on 
more parameters, introducing uncertainty: The estimated optical thickness is 
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influenced by the dust albedo and by the - poorly constrained - extinction 
efficiency, while the thermal radiation relies on the dust temperature. 

Using an extinction efficiency of gext = 1 we find for the optical thickness 
at perihelion at the subsolar point values of r = 0.013 in the insolation 
driven model and r = 0.09 in the jet-model. In the insolation driven model 
the optical thickness of the coma along the line of sight from a point at the 
surface to the Sun is fairly constant over the comet surface and the results 
can well be approximated by the simplified model in [85]. In the jet model 
the optical thickness decreases towards the terminator by a factor of 5. 

Since the extinction efficiency gext is little constrained and can assume 
values in the range gext ~ 1 — 2 we estimate that at perihelion the optical 
thickness at the comet surface is typically in the range r k, 0.01 — 0.03. 
Within jets it can reach values up to r « 0.2. Since the optical thickness is 
proportional to Afp^ the results can easily be scaled to scenarios with different 
levels of cometary activity. 

The intensity received by an observer at 4 km (2 comet radii) from the 
nucleus centre over the subsolar point for different lines of sight is shown 
in Figure [TT] For the insolation driven model, the observed visible intensity 
reaches maxima for lines of sight directed towards the Sun and the nucleus 
(phase angle 180° and 0°, respectively). These peaks are not found in the 
infrared intensity received in the insolation driven model, which shows that 
they are due to the forward and backward scattering peaks of the dust phase 
function (Figure [3]). In the insolation driven model the intensities, both for 
visible and infrared radiation, reach maxima for the lines of sight that touch 
the nucleus tangentially (at about 150°). In the jet-model these peaks are 
not present because most of the dust is concentrated to a narrow region over 
the subsolar point. This also explains the maxima of the infrared intensity 
in the jet-model for the lines of sight in the Sun- and the nucleus-direction 
that are not observed in the insolation driven model. Correspondingly, the 
maxima for the visible radiation on these lines of sight are more pronounced 
in the jet-model than in the insolation driven model. 

The total intensity received by a surface depends on its orientation. To 
give an example, we consider a plane surface at 4 km over the subsolar point 
(see Figure [12]) taking into account also the radiation from the nucleus. The 
intensities emitted or scattered by the nucleus dominate over those from the 
dust; and, since the nucleus is a dark object, the infrared radiation domi- 
nates the total intensity. Note that the nucleus temperature in the model is 
calculated as that of a pure icy surface. Thus the given results for the in- 
frared radiation received from the nucleus are lower limits. If the major part 
of the comet surface is inactive, the surface temperature and consequently 
the infrared flux will be much higher. 
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Fig. 11 Scattered sunlight and thermal infrared radiation at perihelion received 
from a line of sight by an observer above the subsolar point at 4 km from the nucleus 
centre. Only the contribution by the dust is considered, in contrast to radiation from 
the nucleus or direct sunlight. The nominal size distribution a = —3.7 was used. Left: 
insolation driven model, right: jet model. 
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Fig. 12 Radiation at perihelion from all directions on a surface above the subsolar 
point at 4 km from the nucleus centre. Both curves with and without the contribution 
from the nucleus are shown, but no direct sunlight is considered. The nominal size 
distribution a = —3.7 was used. Left: insolation driven model, right: jet model. 



5 Summary and Conclusions 

We have summarised the present knowledge on the dust emitted by Rosetta 
target comet 67P/Churyumov-Gerasimenko and discussed methods and re- 
sults of modching it. The observational data can be described as follows: 
The dust and gas production rates peaked about 30 days after perihelion 
passage during three past apparitions. The coma was characterised by az- 
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imuthal brightness variations probably resulting from active areas on the 
surface. Close to the comet orbit, a distinct line-shaped feature prevailed at 
least from shortly after perihelion 2002 until 2006, beyond aphelion. It is 
interpreted as a pronounced antitail due to the low inclination of the comet 
orbital plane towards the ecliptic. The colour temperatures in both the coma 
and the trail were higher than the blackbody equilibrium temperature at the 
concerned heliocentric distance. 

We have reviewed different approaches to derive parameters of the cometary 
dust production by modelling images of the tail and/or trail, and the results 
obtained by their application to observations of 67P/C-G. We found that 
due to the large diversity of the derived parameter values it is at present 
not possible to formulate a consistent picture of the 67P/C-G dust activity 
and its time evolution. We identified the need to find a common model that 
satisfies all available observations. 

Using the ESA Cometary Dust Environment Model we have predicted 
the terminal speeds of dust, the dust flux in the coma and along sample 
trajectories of the spacecraft, and the radiation flux in the coma. We listed 
results for three mass distributions and five heliocentric distances. We also 
considered different surface activity distributions, among which one implying 
a homogeneous surface of the nucleus and one assuming all dust to be emitted 
by a single active area. The former yielded estimates of quantities averaged 
over space and time, while the latter served to derive upper limits for local 
quantities. 
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